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Abstract 

We investigate the behavior of light baryon masses in three inequivalent large- N c limits: 't Hooft, 
QCDas and Corrigan-Ramond. Our framework is a constituent quark model with relativistic-type 
kinetic energy, stringlike confinement and one-gluon-exchange term, thus leading to well-defined 
results even for massless quarks. We analytically prove that the light baryon masses scale as N c , 
N% and 1 in the 't Hooft, QCDas and Corrigan-Ramond limits respectively. Those results confirm 
previous ones obtained by using either diagrammatic methods or constituent approaches, mostly 
valid for heavy quarks. 
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I. INTRODUCTION 



The study of large- N c QCD finds its origin in the lack of obvious expansion parameter 
for QCD: Setting the number of colors equal to a large value N c instead of 3 allows indeed 
to make perturbative expansions in 1/N C [1]. The relevance of such a framework demands 
that a world in which the number of colors is very large should not be too different from 
our QCD world. The many successes of large-iV c methods in understanding experimental 
as well as theoretical features of QCD provide an excellent a posteriori justification of that 
statement ^j]. Moreover, lattice and AdS/CFT-based calculations of various observables 
using different gauge groups strongly favor the idea that SU(3) is actually close to SU(oo) 
(see e.g. Q). 

The standard way of generalizing QCD to a large number of colors is also the first one 
proposed in the literature, by 't Hooft It consists in taking the quarks in the fundamental 
representation of SU(iV c ), and then letting the number of colors becoming arbitrarily large 
while the so-called 't Hooft coupling g 2 N c remains constant for any iV c (g is the strong 
coupling constant). Note that in the 't Hooft limit, the number of quark flavors remains 
finite; later Veneziano would propose to set Nf = 0(N C ) Q] so that the theory remains 
planar, but the internal quark loops are no longer suppressed as in the 't Hooft limit. 
However, it has been noticed soon after 't Hooft proposal that the extrapolation of QCD 
to arbitrary N c is not unique: The only theoretical constraint is indeed that the resulting 
theory reduces to QCD at N c = 3 In this last reference for example, it has been suggested 
that some quark flavors could be in the ^4at c -2 representation, that is the (N c — 2)-indices 
antisymmetric representation, equal to the fundamental one at N c = 3. This is called 
the Corrigan-Ramond limit. Moreover, since the fundamental and A 2 representations are 
equivalent for SU(3), another large- N c limit can be proposed in which quarks are taken to 
be in the A 2 representation of SU(iV c ). Denoted QCDas, that limit interestingly leads to a 
theory equivalent to M = 1 SUSY Yang-Mills, as shown in 

In any large- N c limit, mesons and glueballs are always quark- ant iquark and gluonic bound 
states respectively, but baryons need to be considered more carefully. As pointed out by 
Witten {7 1, baryons in 't Hooft limit should be seen as bound states made of N c quarks 
in a totally antisymmetric color singlet. In that reference, using Feynman diagrams and a 
nonrelativistic quark model, it has been shown that the baryon masses scale as N c at the 
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dominant order. It was also suggested that the same result should hold for baryons made 
of light quarks. Remark that corrections in 1/N C can be added to improve the agreement 
of baryonic mass formulas with experiment (see e.g. |8j). In the Corrigan-Ramond limit 
however, a baryon can be made of three quarks as in QCD, so its mass should stay constant 
with N c at dominant order j^]. Finally, baryons in the QCDas should rather consist in 
N C (N C — l)/2 quarks in a totally antisymmetric color singlet [9], with a mass scaling as N% 
as shown in [lo| by using diagrammatic methods similar to the ones of 

In the present paper, we propose a constituent quark model that has the peculiarity of 
being relevant in the light baryon sector: Quarks can be massless thanks to a kinetic term 
of relativistic form, and long-range confining interactions are taken into account within a 
flux-tube (or string) picture. Our model is discussed in Sec. Ull where we keep the formalism 
general so that any baryonic system can be studied, either in QCD or in the aforementioned 
large- N c limits. Then, using mathematical tools that are detailed in Appendix |AJ we find 
analytical upper and lower bounds for the light baryon mass spectra in the 't Hooft, QCDas 
and Corrigan-Ramond limits in Sees. Illlj. IIVI and |V] respectively. As finally summarized in 
Sec. I VI} those bounds allow to find how the light baryon masses scale with N c and provide 
a confirmation of results obtained in previous works devoted to that topic. 



II. BARYONIC HAMILTON! ANS 

A. General considerations 



A Hamiltonian describing bound states of baryonic type is needed as a starting point for 
our study. Since the gauge group considered is SU(iV c ), and since the quarks are allowed 
to be in an arbitrary color representation R of that group, we mean by "baryonic type" a 
bound state made of quarks only, whose color function is a totally antisymmetric singlet 
with respect to the exchange of two quarks. Let us denote n q the minimal number of quarks 
necessary to build such a color singlet. Then the dominant term in the quark kinetic energy 
can be written under the spinless Salpeter form 

n q 

i=i 
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which has the advantage of being well-defined even for massless quarks. It is sometimes 
called semirelativistic in the literature since it is not a covariant formulation. 

Let us now discuss the interactions between the quarks within a generic baryonic system. 
At the perturbative level, i.e. at short distances, one-gluon-exchange processes are dominant. 
The corresponding short-range potential can be computed from the QCD Feynman diagrams 
at tree level. In the quark-quark case, one is led to the well-known Fermi-Breit interaction, 
whose leading part is spin-independent and of Coulomb form. Typically, one expects the 
one-gluon-exchange term to be 

Voge = \{C n - 2C q ) a s jr _L_, (2) 

I 1. *" J I 



i<j=l 

,2 



where a s = g /4ir and where C qq and C q are the quadratic Casimir operators of SU(iV c ) in the 
representations of the quark-quark pair and of the quark respectively. Xi obviously denotes 
the position of quark i. Notice that we assumed that each quark pair in the system is in the 
same color channel. The reason of such a choice will appear more clearly in the following 
sections. As previously mentioned, the 't Hooft large- N c limit is such that N c — y oo and 
a s = ao/N c , with ao independent of N c . Since a s ~ 0.35-0.40 in various potential models of 
real QCD, q;q ~ 1. With such a value, all arguments within the square roots appearing in 
the following to compute baryon masses are well definite positive numbers. 

The long-range part of the interactions, responsible for the confinement, will be described 
within the framework of the flux-tube model. From the Casimir scaling hypothesis, one 
considers that each color source generates a straight flux tube whose energy density (or 
tension) is proportional to its quadratic color Casimir operator. Then, the flux tubes have 
to meet in one or several points such that the total energy contained in those flux tubes 
is minimal. In "QCD" baryons (three quark systems), the junction point is identified with 
the Steiner (or Fermat or Toricelli) point of the triangle made by the quarks, leading to a 



Y-junction configuration as observed in lattice computations llMl3|. Excepted for highly 



asymmetric quark configurations, that are not expected to be relevant for low-lying baryons, 
the junction point can be identified to the center of mass of the system within an excellent 
approximation [l^j]. Inspired by the usual baryonic case, we will assume a confining potential 
of the form 

V c =^aJ2k~R , (3) 

° n i=i 1 
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where R is the center of mass of the system. According to the Casimir scaling, (C q /C D )a is 
the flux-tube energy density since C D is the quadratic Casimir in the fundamental represen- 
tation and a is the fundamental string tension, that is the one between a quark-antiquark 



pair. It is expected that a is constant with N c at leading order 15) . In various potential 
models, a ~ 0.15-0.20 GeV 2 . Actually, there may be more than one junction point for 
multiquark states (tetraquarks, pentaquarks, ...), and algorithms exist to compute them 
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10] as well as lattice computations |17|). So, a baryon made of //,, quarks is probably 



a very complicated spatial structure with numerous strings connecting quarks and multiple 
junctions in order to minimize the potential energy. Nevertheless, besides the fact that 
potential (EJ) is of one-body form, which is better in view of analytical calculations, it has 
several interesting physical features: 

• When N c = 3, it reduces to an excellent approximation of the genuine Y-j unction. 

• Provided that a single junction point is present, its motion can be neglected at large- 



N c (see arguments in [?]) and its position can thus be identified with the center of 
mass. 

At large- N c , the string tension only depends on the iV-ality of the representation (in our 
case, the number of fundamental representations needed to build a given representation 
by tensor product). This has been shown using several approaches such as holographic 



techniques or lattice calculations [18l . Il9| . If k fundamental strings emerging from k 
quarks in the fundamental representation connect to a junction, the energy density 
of the string emerging from this junction is expected to be = k a at the dominant 
order in N c . The corresponding energetic cost is similar to the one for k fundamental 
strings. So, the connection of the strings emerging from all quarks in the baryon 
into a single point is not energetically disadvantaged with respect to the apparition of 
multiple junctions inside the baryon at large- N c . 

For these reasons, we will keep potential ([3]) in our calculations, since it seems to be relevant 
for both iV c = 3 and N c — > oo. 

It is worth mentioning that the total Hamiltonian, defined by 

H b = T + V oge + V e , (4) 
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is very similar to the string model of light baryons at large N c proposed by Witten in his 
seminal paper Q. The peculiarity of the present work is that the model will be generalized 
to other large- N c limits than 't Hooft's one and that we use techniques allowing to obtain 
approximate analytical mass spectra for Hb- A last comment that can be done about our 
model concerns the strong coupling constant. Already at one-loop, as is running with 
the energy scale q 2 and we know that a s (q 2 — > oo) = 0. Moreover, we consider that 
a s{<i 2 0) tends to a finite, nonzero value that corresponds to a s used in This is 
coherent with the aforementioned lattice QCD calculations, showing in general that the 
static potential between quarks can be accurately fitted by a potential separated into a 
dominant nonperturbative part (the flux tubes) and a residual perturbative part (the one- 
gluon-exchange potential with a "frozen" value of a s ). 

The spin effects will be completely ignored in this work. For instance in the real world, 
the ratio (jtoa — fn N ) /\{m^ + m N ), which is attributed to spin (possibly isospin) dependent 
interactions, is 27% (around 1/N C with N c = 3). This spin contribution is not negligible 
but can be considered as a perturbation of order 1/N C in the large- N c limit js]. Previous 
works have considered the 1/N C expansion for the mass operator in terms of various spin or 



isospin dependent operators inspired by simple constituent quark models 
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2l|. Explicit 



calculations of the spin contribution have been calculated for N c = 3 in the framework of 



a nonrelativistic quark model 



21] . This problem for a relativistic kinematics and arbitrary 



values of N c can be treated within our formalism and will be studied in a subsequent paper. 
Here, we only focus on the dominant effects in baryons. 



B. Upper and lower bounds 



Summarizing the discussion of the previous section, we are left with the baryonic Hamil- 
tonian 



i=l 



X; — R 



E 



t<j=l J 



where a = ^-a, and b = ^(2C 9 - C qq )^-. 



(5) 



(6) 



Since we focus on light baryons, we have set rrii = 0, which is a good approximation for the 
u and d quark masses. 

Although finding the exact mass spectrum of Hb, denoted M&, would require numerical 



6 



computations, analytical techniques allow to find upper and lower bounds of that spectrum. 
For the sake of clarity, we leave the technical details for Appendix |A] while we give here the 
final results. From (IA2j) and (IA12j) . we can deduce that the mass Mb of a given eigenstate 
is such that 

1 



n„ inf ( ib 

m \ 



v + —, 



a b 
n„ r 



V>) <M b < Jia(Qn q -bP 3 / 2 ), 



(7) 



where we recall that P is the number of quark pairs and Q is the band number of the consid- 
ered state in a harmonic oscillator picture. The upper bound, obtained by the auxiliary field 
method (AFM), is state dependent, while the lower bound only concerns the ground state 
(hence, it is also a lower bound of the whole spectrum). Following flA13j) . it can be approx- 
imated by a/2glP (2 — b{n q — 1)). The number of excitation quanta, K = Ylili {^ n i + Pi), 
could be seen as either of order 1 or of order n q . We will consider in the following that 
K = 0(1) in order to get simpler mass formulas, but we point out that our final results 
would remain valid if K = 0(n q ). 

Whatever the values of a and b, ([7]) and (1A3|) implies that 



Ml <aK + /3, 



(8) 



where a and (3 do not depend on quantum numbers. So, for a fixed value of n q , the upper 
bound predicts a Regge behavior for orbital and radial excitations. This is observed in the 
real QCD world with N c = 3. Moreover, large- N c methods help to understand why baryonic 
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and mesonic Regge slopes are equal 

A mean field approximation has been used in various papers [?,l9| to estimate the baryon 
masses in the large- iV c limit. One can ask what are the connections of our model with this 
approach? The energy e of a single quark in a central mean field with a funnel form ar — jH/r 
is given by 

e = v/4a(Q* - 0) with Q = 2n* + I* + 3/2, (9) 
within the AFM approximation 



23]. The mass M* of n q independent quarks is 



m* = J2 VMQt - Fj. 



(10) 



For large values of n q , we can assume a good equipartition of the energy Q* Q/n q , where 
Q is the number of quanta for the baryon state considered. The corresponding mass is then 



M* ^ J4a(Qn q - 0nf, 



(11) 



With the identification a = a and (3 = 2~ 3 ^ 2 bn q , the upper bound ([7]) is recovered. Never- 
theless, a mean field approximation is not fully equivalent to our approach. For instance, 
in the mean field formulation, the wavefunction is the product of n q individual functions 
4>j(xj — R), which is very different from the wavefunction flA6p . 



III. 'T HOOFT LIMIT 



In this case, the quarks are in the fundamental representation of SU(iV c ), and the number 
of colors becomes large while the 't Hooft coupling, or equivalently a s N c , remains constant 
for any iV c . The number of quark flavors is finite. The generalization of a baryon is then 
given by a state with n q = N c , the N c quarks forming a totally antisymmetric color singlet, 
implying that any quark pair is in the A2 representation B jz|- Thus the Casimir operators 
read in this case 



C„ — Cr 



N 2 c -l 

2N r . 



Cqq - 



(N C -2)(N C + 1) 



(12) 



The interested reader can find in e.g. 



24J a way to compute the Casimir operators of SU(iV c ). 



This implies that the a and b factors fl6]) are now given by 



a = a, 



N c + 1 
2N? 



-a . 



(13) 



Due to the large value of JV C , one has Q -> 3N c /2 (K = 0(1)), P ->■ iV c 2 /2, and © becomes 



N c inf ( ib 



'p* 



2 r ~4r 



r)< } fi, < .Y, \ / <t ( 6 — 



(14) 



Following flA13j) . this last lower bound is close to the explicit formula N c y a (2 — More- 
over, if n s strange quarks with nonzero mass are present, it is readily obtained from (1A4j) 
that they bring a mass term approximately equal to 

,2 



AM, w n 



(15) 



Using ( 1A8j) . the mean square radius of the baryon, at the limit N c — > 00, is given by 

/ 2 \ 6 — \/2«o 



x 1 4a 

The size of the baryon tends toward a finite value in this limit, as suggested in 



(16) 
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Since a and «o do not depend on N c , (j!4j) actually shows that any eigenvalue of Hamil- 
tonian (jSj) scales as iV c in the 't Hooft limit. We argued in Sec. [II] that this last Hamiltonian 
provides a relevant description of baryonic systems, at least at the dominant order; conse- 
quently, we have shown that the light baryon masses scale as 

M b oc N c (17) 



in the 't Hooft limit, just as it has been shown for heavy baryons in [2j . This is in agreement 



with the results in 



0. 



251 ] which show that this scaling must be verified whatever the number of 



26J . Moreover the mean 



light quark flavors. The strange quark contribution (|15|) is constant with N c and proportional 
to m^, which is the SUj?(3)-breaking parameter in our model. That behavior is indeed 
coherent with former results obtained in large- N c baryon physics 
square radius is also of order 0(1). 

We point out that, although Q (or K) does not appear explicitly in the formulas anymore, 
these results are valid for any excited state of Hj, as well as for its ground state. It can finally 
be observed, from the formulas of Appendix |AJ that a pairwise confining potential would 
not have led to masses of order N c excepted if a was of order 1/N C , which seems not relevant 
as mentioned above. 



IV. QCDas LIMIT 



That limit is different from the 't Hooft one since the quarks are now assumed to be 
in the A2 representation B of SU(iV c ), demanding the baryonic color wave functions to be 



ly antisymmetric color 
9j. As for the 't Hooft 



reconsidered. Indeed, the number of quarks needed to build a tota 
singlet at any iV c is now n q = N C (N C — l)/2 as explicitly shown in 
limit, n q is equal to the dimension of the quark color representation, meaning that all the 
quark color configurations appear only once in the singlet. If it is readily computed that 
C q = Cg = (N c — 2)(N C + 1)/N C , it is nevertheless worth discussing a bit the value of C qq . 
A quark pair inside the baryon can be in the following color channels: 



(18) 



However, the last channel is the only relevant one, for the following reasons. From a quark 
model point of view, one can check that EE3 leads to repulsive interactions, that are thus not 



appropriate to describe a baryon. Moreover, the totally antisymmetric channel is forbidden 
from a diagrammatic point of view, as shown in j^. The point is that a quark now carries 
two fundamental color indices. Consider a one-gluon-exchange process involving initially 
two quarks with indices [ab] and [cd] . The four indices must be different in order for the pair 
to be in the A4 representation. The exchanged gluon would then change the pair into, say, 
[ad] and [bc\. But, in virtue of the Pauli principle, that pair is already present in the baryon, 
leading to the vanishing of that process. From that argument, the conclusion is that any 
quark pair must have only one common color index, typically [ab] and [ac] . This corresponds 
to the last diagram of (fl8"l) . and one finally concludes that C qq = Cm = 2(N% — N c — 4)/N c . 
The a and b factors are given, in the QCDas case, by 

2(N c -2) ,2 

a = iv c -i a ' iVf" ' (19) 

and, as N c ->■ 00, it can be checked that Q -> 3iV c 2 /4 (K = 0(1)) and P ->■ iV 6 4 /8. 
Equation (j7j) becomes 



2 

c 

2 



N' 

c 



inf (*\y& + (<,r- I V) < M b < N 2 Ja (3 - . (20) 



Following (IA13[) . that last lower bound is approximately equal to A^^ya (l — 9 f). The 
strange quark contribution is moreover given by 



,2 



while the baryon radius reads 

CO - (22) 

Using the same arguments as in the 't Hooft limit, we can conclude from (120]) that we 
have shown that the light baryon masses scale as 

M b oc N 2 C (23) 

in the QCDas limit, as proposed in 28] and proved in Witten's diagrammatic way in flo| |. 
Moreover, the SUi?(3)-breaking term brings a contribution independent of N c in this case 
also, as pointed out in j^. The behavior is similar for the baryon radius. Notice that, 
with quarks in the A 2 representation, a totally symmetric color singlet could be built with 
iV c quarks, that would lead to a baryon mass scaling in following This can be 
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understood within our framework also, because the spatial wave function should then be 
antisymmetric to preserve Pauli's pri ncip le. In that case, it can be computed within our 



framework that Q oc n^ 3 at large n q 29j, leading to Mf, oc Nc /0 following (IA2j) and (r 2 ) oc 



7/6 



Nc /S following (IA8|) . 

V. C O RRI G AN- RAM O ND LIMIT 

In that limit, baryons are three-quark states for any N c : Two quarks are in the fundamen- 
tal representation and one (say quark number 3) is in the Aw c -2 one. The total interaction 
potential of our model reads, in the limit where N c — > oo and where the 't Hooft coupling 
is fixed, 

+ V c a - R\ + |fi -R\ + \x 2 -R\)-^ ( F 1 _ . + F 1 - ■ ) . (24) 

\ / Z \\Xi — £3 1 p2 — ^3|/ 

That potential, being constant with iV c , would lead to constant masses also since n q = 3. 
However, it is rather problematic since the baryon picture suggested is quite far from what 
is known at N c = 3. In particular, no one-gluon-exchange term is present between the 
quark 2 and 3: In the previous limits, the 1/N C scaling of this term was compensated by the 
number of quarks scaling as iV c at least, but here n q is finite. Moreover, the asymmetry of 
the confining potential would lead to a tower of excited states quite different from what is 
expected. For those reasons, we prefer not to investigate more the Corrigan-Ramond limit. 
We mention for completeness that the phe nomenology of the Corrigan-Ramond limit has 



already been studied in previous works [28|, |30(, and found to be less realistic that the one 



coming from the 't Hooft or QCD as limits, in agreement with our discussion. Moreover 
another limit which is somewhat in between the 't Hooft and Corrigan-Ramond ones has 
also been proposed {31], but it is out of the scope of the present work since it requires a 
formalism in which quarks are Dirac spinors. 

VI. CONCLUSIONS 

In this paper, we have built a constituent model describing a baryonic system, that is 
a bound state made of quarks only. Our approach is a generalization of Witten's original 
proposal |Z|, i.e. a one-gluon-exchange potential plus a stringlike confinement, but in which 
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a kinetic term of relativistic form is added. Using analytical methods, we have been able to 
find analytical upper and (approximate) lower bounds for the baryon mass spectrum in the 
limit where quarks are massless. This is a peculiarity of the present work since analytical 
results are generally known only for heavy quarks. 

The most important point is that we have confirmed that Mb = 0(N C ) for light baryons 
in the 't Hooft limit, and that M& = O(N^) in the QCDas limit in agreement with the recent 



work 



10j. Moreover, M& = 0(1) in the Corrigan-Ramond limit, but this last case does not 



seem relevant for a constituent approach. We have checked that these results are obtained for 



the more general confining potential a 



Xj — R 



with p > 0. But, whatever the value 



of N c , baryon masses lie on Regge trajectories, for orbital and radial quantum numbers, only 
for p = 1 . 

The contribution coming from strange quarks is of order 0(1) in the three considered 



cases, as suggested in [281 ] . It can also be checked that the typical baryon size is independent 
of N c for the three considered limits, recovering a result obtained within a Skyrme model of 
baryons jjj]. 

Those results can basically be summarized as follows: For light baryons, M& is propor- 
tional at the dominant order to the number of quarks needed to build a totally antisymmetric 
color singlet. This seems rather intuitive a posteriori but it was actually not trivial a priori 
for the following reasons. First, the quarks are massless, and the usual guess Mj, oc n q m q is 
not applicable. Second, the explicit form of the potential does not allow to deduce such a 
scaling law. Third, a simple counterexample can be found as shown in Sec. IIV1 If the color 
singlet is totally symmetric, the baryon masses do not scale as n q at large n q , showing the 
nontrivial character of our result. 
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Appendix A: Useful formulas 



In this appendix, we give analytical approximate mass spectra for a generic Hamiltonian 
of the form 



i=l 



pf + mf + at 



Xi — R 




n q 


a 2 \x t 


x j 




i<j=i 







(Al) 



in view of an application to baryonic systems studied in the present paper. First of all, 
we focus on the ultrarelativistic limit rrii = 0. Although standard quantum mechanical 
techniques have difficulties to handle that limit, the auxiliary field method (AFM) has been 
shown to be efficient in that case also. We refer the reader to 2^| for a detailed discussion 
of that method, but for our purpose it is sufficient to mention that it leads to the following 
upper bound for the mass spectrum of Hamiltonian (1A1[) 



M 2 ul {Q) = 4( 0l + a 2 P 1 / 2 )(Qn g - 6P 3 / 2 ), 



with 



n a — 1 



Q 



-(n a 
2 V q 



l)+K, K=J2 i 2 ^ + 4)> and p 



n q (n q - 1) 



(A2) 



(A3) 



i=l 



Q is the band number (in a harmonic oscillator picture) for the eigenstates of H, and P is 
the number of particle pairs. It is assumed that the value of b is such that is a positive 
number. In the case where a\ = or a 2 = 0, (IA2I) reduces to formulas previously found in 



29| | . Notice that, if n s (< n q ) quarks have a small nonzero mass m s , the corresponding mass 
shift is given by 2^| 

2 

771 

AM s (Q)*n s —f- (A4) 
where fMo(Q) is a massless quark's average kinetic energy, given by 



QM ul {Q) 



(A5) 



Qn q -bP 3 / 2 2(Qn q -bP*/ 2 y 
The AFM not only leads to approximate analytical mass spectra, but also to approximate 



analytical wave functions, a 



Combining results from 



29 



lowing to compute various observables with a good accuracy [32 ]. 



32] , a n 9 -body eigenstate is written 

n q — 1 



(A6) 
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where ip n -£-{\j,yj) is a three dimensional harmonic oscillator wave function, depending on 
the Jacobi coordinate y*j 29] and decreasing asymptotically like e~ A ?^ 2 . The quantum 
numbers are such that Y^j=i (^ n j + ^j) = K an d the scale parameters Xj are given by 



Xn 



In 



j + 1 V Q 



(A7) 



The state (IA6I) has neither a defined total angular momentum nor a good symmetry, but its 
is characterized by a parity (— 1) . By combining states ( 1A6I) with the same value of K (or 
Q), it is generally possible to build a physical state with good quantum numbers and good 
symmetry properties, but the task can be technically very complicated [33 ]. 



The baryon mean square radius can 
to the mass spectrum M nl , leading to 



je computed with the eigenstates (IA6[) corresponding 



29 



32 1 



/ 1 nq 



Rf 



Q 



n q /j, (Q) 



(A8) 



Since this result only depends on the quantum numbers via Q, and since a physical state 
must be a combination of eigenstates with the same value of Q, (1A8|) is also valid for a 
physical state. 

In order to check our calculations, an upper bound for the ground state of Hamilto- 
nian (lAlj) has been obtained by a variational method with a trial function given by ( 1A6|) 
with K — and new scale parameters Xj = \[JJ 
minimizing (<f)\H \ <p) with respect to A (see e.g. 



7 + T)A, where A is a free quantity. By 



29 



34] for a computation of the needed 



matrix elements), this leads to the ground state upper bound 

16 

M u2 = —P{ai + a 2 VP){2 - bVP). 

7T 



(A9) 



Not only upper bounds can be obtained, but also a lower bound. Remarking that 

n q n q 

y y j \xj - r > — ^2 %j\ ( Ai °) 



| *^ J | 

i=l " q i<j=l 

thanks to the triangular inequality, eigenenergies of the Hamiltonian 

b 



i=l i<j=l L 



X j Ob <j 



(All) 
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with a = a 2 +a 1 /n q , are lower bounds of the eigenenergies of Hamiltonian (lAip . Furthermore, 
a lower bound on the ground state mass of H 1 can be found in |35|, that is 



Mr 



n, 



inf ( ib 

\4>) \ 



a/p 2 + m 2 + 



ar 



4 



(A12) 



More explicitly, an approximate formula for Mn at m = can be found using the AFM 



23 1 



Mji » M? 2 = 2P(a 1 + a 2 n q ) [2 - b(n q - 1)] 



(A13) 



The variational character of M\ 2 cannot be guaranteed, but since this last AFM approxima- 
tion is very accurate |23|, it can be reasonably supposed that M i2 is still a lower bound of 



the ground state mass of (IA1I) . This has been numerically checked on several systems for 



n q = 2 and 3. 



Since min(Qn q ) = 3P, one can verify that M U \(Q) > M u2 > M\ 2 for n q > 2, as expected. 
The crucial point of these calculations is that the three bounds have the same behavior for 
large values of n q . So we can have confidence on the large n^-behavior of the exact solutions. 
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